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The Dicke model describes an ensemble of N identical two-level atoms (qubits) coupled to 
a single quantized mode of a bosonic field. The fermion Dicke model should be obtained by 
changing the atomic pseudo-spin operators by a linear combination of Fermi operators. The 
i— i . generalized fermion Dicke model is defined introducing different coupling constants between the 
>■ ; single mode of the bosonic field and the reservoir, gi and g 2 for rotating and counter-rotating 
terms respectively. In the limit iV — > oo, the thermodynamic of the fermion Dicke model can 
be analyzed using the path integral approach with functional method. The system exhibits a 
second order phase transition from normal to superradiance at some critical temperature with 
the presence of a condensate. We evaluate the critical transition temperature and present the 
spectrum of the collective bosonic excitations for the general case (gi ^ and g2 ^ 0). There is 

^ \ quantum critical behavior when the coupling constants g\ and g 2 satisfy g\ + gi = (u u) a, where 
u is the frequency of the mode of the field and Q is the energy gap between energy eigenstates of 
the qubits. Two particular situations are analyzed. First, we present the spectrum of the collective 
bosonic excitations, in the case g± ^ and #2 = 0, recovering the well known results. Second, the 
case g\ = and gi ^ is studied. In this last case, it is possible to have a superradiant phase 
when only virtual processes are introduced in the interaction Hamiltonian. Here also appears a 

quantum phase transition at the critical coupling g 2 = (wo Q) a , and for larger values for the critical 
coupling, the system enter in this superradiant phase with a Goldstone mode. 
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1 Introduction 



The aim of the present paper is to investigate the thermodynamic in a generalization of the 
Dicke model [1], using the path integral approach with the functional integration method. We are 
considering the question of how does the counter-rotating terms of the interaction Hamiltonian 
contribute to the system exhibits a phase transition from normal to superradiance at some critical 
temperature f3~ l . 

We are particularly interested in consider N identical two-level atoms interacting with a quan- 
tized bosonic field, i.e., the spin-boson problem. The most discussed situation in the literature is a 
system of atoms interacting with a reservoir which consists of a surrounding quantized electromag- 
netic field. In this situation the frequencies associated with the electromagnetic field are continu- 
ously distributed and the atomic systems are characterized by a denumerable set of frequencies. 
Nevertheless, different situations has been extensively studied. One is the case of electromagnetic 
field in cavities where the field modes are described by a set of denumerable frequencies [2] [3]. 
One particular example is the Jaynes-Cummings model [4], where was investigated the dynamics 
of a single two-level atom interacting with one mode of a quantized bosonic field in an ideal lossless 
cavity. The eigenstates of the model, i.e., the dressed states of the two- level atom-bosonic mode 
were first obtained by these authors. After these important papers, the spin-boson problem has 
been studied by many authors, receiving considerable discussions. The development of new tech- 
niques has lead to increase the interest in spontaneous radiations effects in such systems. More 
recently, this kind of model has been used also to analyze dissipation in quantum computers. A 
model describing a two-level atom coupled to a reservoir of harmonic oscillators was introduced 
by Di Vicenzo [5], to study decoherence in quantum computers. This model presents destruction 
of quantum coherence without decay of populations, owing to a quite particular coupling between 
the qubit and the reservoir. 

In the following we are discussing the Dicke model, where a single quantized mode of a bosonic 
field interacts with a reservoir of N identical two- level atoms (qubits) at temperature One 
method to obtain information from the small system of interest, eliminating the reservoir degrees 
of freedom is to use the master equations, introducing the reduced density operator [6] [7]. An 
alternative and non-perturbative powerful tool to deal with the elimination of the reservoir degrees 
of freedom, is the path integral approach with functional integration method [8]. In a general 
situation, this method allow us to obtain the semiclassical equation of motion and also gives a 
systematic way for the evaluation of quantum corrections. 

Although functional integrals have been used in different area of physics, as for example quan- 
tum field theory and also statistical mechanics, in the past this technique have no found much 
use in quantum optics. An exception is the Zardecki paper [9]. Other paper that deserves to 
call attention is the Hillery and Zubaire [10] work. In this paper the path integral approach in 
quantum optics was developed, where systems with a few bosonic modes were studied. Using 
the Schrodinger-Glauber-Sudarshan's coherent states [11] [12] [13] [14] [15] [16] representation of 
these modes, these authors studied the propagator in a single mode system and also the case of N 
modes of a bosonic field. Since, in the Dicke model, the coupling between the N qubits and the 
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single quantized mode of the bosonic field is linear, in order to analyze the nonanalytic behavior 
of thermodynamic quantities, the integration over the degrees of freedom of the environment can 
be performed exactly. 

To apply the path integral approach with functional method, as a technique to investigate the 
thermodynamic of the Dicke model, two steps are mandatory. First, it is necessary to change the 
atomic pseudo-spin operators of the model by a linear combination of Fermi operators to define 
the fermion Dicke model [17] [18]. A generalization of this model can be achieved introducing 
the counter-rotating terms in the interaction Hamiltonian. Second, the thermodynamic limit 
(N — > oo), where N is the number of qubits must be taken [19] [20]. Therefore, studying the 
fermion Dicke model with counter-rotating terms, describing virtual processes, we are interested 
to evaluate the critical transition temperature of the model. In this generalized model, one may 
ask the importance of the contribution of each of these processes, real ones and virtual ones in the 
formation of the condensate. 

A point that it is worth commenting upon is the relationship between our questions concerning 
the contributions of the real and virtual processes in the superradiant phase transition and how 
does the vacuum fluctuations and the self-reaction contribute to the radiative processes in atoms. 
Many authors [21] [22] [23] used the Heisenberg picture to identify the vacuum and the self-reaction 
contributions in radiative processes in atoms. Although there was some attempts to separate the 
contribution of the vacuum fluctuations and radiative reaction in the spontaneous decay [24] [25], 
Dalibard et al [26] [27] discussed how does the magnitude of these separate effects can varied by 
means of the particular ordering chosen for atomic and creation and annihilation field operators. 
When normal ordering is used, the source field gives the total contribution to the decay, and 
for anti-normal ordering, both the source field and vacuum field gives equal contributions to the 
spontaneous decay. Finally, if the symmetric ordering is chosen, the vacuum and the source field 
gives different contributions for the decay. There is no ordering where the spontaneous emission 
can be entirely associated to the vacuum field [28]. 

As we already emphasized, we are interested to study the nonanalytic behavior of thermody- 
namics quantities near a phase transition in the fermion Dicke model where the counter-rotating 
terms are also presented. Introducing different coupling constants, g\ and Q2 for rotating and 
counter-rotating terms respectively, we are able to identify the contribution of each of the pro- 
cesses, real ones and virtual ones in the formation of the condensate. The order parameter of 
the transition is the expectation value of the number of excitations associated to the mode of the 
bosonic field per atom. We evaluate the critical transition temperature of the model and present 
the spectrum of the colective bosonic excitations in the general case, for the case gi ^ and 
(?2 = and also g 1 = and gi ^ 0. In both cases it appear a critical behavior with Godstone 
(gapless) modes. In the last case, it appears a quantum phase transition [29] at critical coupling 

(?2 = (^o^) 5 - Therefore, in the fermion Dicke model with the counter-rotating terms it is possi- 
ble to have a quantum fase transition and a superradiant phase, where the radiation rate on the 
number N of atoms becomes a quadratic dependence [30] , in opposite to the situation of N atoms 
radiating incoherently, where the radiation rate is proportional on the number of the atoms N. 



2 



At this point we would like to make a summary of results concerning the thermodynamic 
of the Dicke model. An important result were obtained by Hepp and Lieb [31]. These author 
presented the free energy of the model in the thermodynamic limit. For a sufficiently large value 
for the coupling constant between the qubits and the single quantized mode of the bosonic field, 
the model present a second order phase transition from the normal to the superradiant phase. 
Latter, using a coherent state representation they generalize some results for atoms of more than 
two-levels without assume the rotating-wave approximation. They also investigated the stability 
of the model with an infinite number of bosonic modes [32]. The study of the phase transitions 
in the Dicke model was presented also by Wang and Hioe [33], where some of the results obtained 
by Hepp and Lieb were rederived. The generalized Dicke model, where the counter-rotating terms 
are also present in the interaction Hamiltonian, was investigated by Hioe [34] and also Duncan 
[35] . Hioe studied the thermodynamic of the generalized Dicke model with two different coupling 
constants using the Schrodinger-Glauber-Sudarshan's coherent states. Duncan proved that the 
critical temperature of the generalized model is higher than the temperature of the original model. 
Also a bosonization procedure was employed to study the phase transitions in the generalized 
Dicke model. Employing a Holstein-Primakoff mapping [36] [37] [38], which express the angular 
momentum in terms of a single bosonic mode, Emary and Brandes [39] [40] were able to express 
the generalized Dicke model in terms of two mode bosonic field. These authors discussed the 
relation between the quantum phase transition and the chaotic behavior that appear in the model 
for finite N. An indication of quantum chaos is the change in energy level-spacing statistics from 
Poissonian to being described by the Gaussian ensemble of random matrix theory. This chaotic 
behavior was discussed also in the Jaynes-Cummings model by Graham and Hoherbach [41] [42] 
and Lewenkopf et al. [43], in the situation where the counter-rotating terms are present in the 
interaction Hamiltonian [44] [45], since the seminal paper of Milonni et al. [46] 

We would like to stress that we are not doing any distinction between the elementary unit 
of quantum information, the qubit and the two-level atom. Our intention in this paper is to 
review some important spin-boson models and also the basic ideas of the path integral approach 
with functional methods that can be used in these systems. This paper is organized as follows. 
In section II we discuss two-level atoms-boson field interaction Hamiltonians. In section III the 
Hilbert space for the fermion Dicke model is discussed. In the section IV the path integral with 
functional integral method is applied to study the thermodynamic of the generalized fermion Dicke 
model. We evaluate the critical transition temperature of the model for the general case (g\ ^ 
and g 2 ^ 0) and also present the spectrum of the collective bosonic excitations, for the case g 1 ^ 
and g 2 = and also for the case g x = and g 2 ^ 0. Conclusions are given in section V. In the 
paper we use ks — c — h — 1. 
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2 The two-level atoms-Bose field Hamiltonians 



In order to describe the dynamics of the reservoir and the small system, we have to introduce 
the Hamiltonian governing the interaction of the quantized bosonic field with free non-identical 
two- level atoms. Free means that there is no interaction between the two-level atoms (qubits). We 
are assuming that the qubits are enclosed in a very large lossless cavity, and the distance between 
the qubits is large enough, so that the interaction between the qubits can be neglected. 

The main purpose of this section is to review a number of relevant qubit-Bose field Hamiltoni- 
ans. This discussion is standard in the literature. For a recent treatment see for example the Ref. 
[47]. It is worth mentioning that there are a variety of theoretical models of reservoirs. We have 
chosen the models that seems to us most significant and interesting. The first situation is when 
the system S is coupled to an infinite number of harmonic oscillators. In this situation there are 
two kinds of reservoir of common interest. The first one is a thermal reservoir, where we assume 
that the harmonic oscillators are in thermal equilibrium at temperature f3~ l . The second one is 
a squeezed reservoir. The specific system-reservoir model which is appropriate for the study of 
several interesting situations is when the harmonic oscillator bath is constituted by a bosonic field 
in free space also in the presence of macroscopic structures. 

Therefore let us consider a bosonic quantum system S, with Hilbert space which is coupled 
with the reservoir of qubits, with Hilbert space 1iS B \ Let us assume that the reservoir is in thermal 
equilibrium at temperature The bosonic quantum system is a sub-system of the total system 
living in the tensor product space 

Let us denote by Hs the Hamiltonian of the quantized bosonic field (where the "S" means 
small system), by Hb the free Hamiltonian of the iV-qubits (where the "B" means bath) and Hi 
the Hamiltonian describing the interaction between the quantized bosonic field and the qubits 
reservoir. The Hamiltonian for the total system can be written as 

H = H S <g> Ib + Is ® H B + Hi, (1) 

where 1$ and 1b denotes the identities in the Hilbert spaces of the quantized bosonic field and the 
qubit reservoir. 

The free j — th qubit Hamiltonian will be denoted by Hp\ Therefore we have 

tf^lOi^lOi, (2) 

where | % )j are orthogonal energy eigenstates accessible to the j — th qubit and uo\^ are the 
respective eigenfrequencies. Using Eq. (2) and the orthonormality of the energy eigenstates we 

can write the j — th qubit Hamiltonian Hp* as 

^=j:^(\i){i\) . (3) 

i=i v J o 
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Throughout the article, for the sake of simplicity in the notation, we are assuming that the operator 
= and the same convention is used for the other operators. 

Let us define the pseudo-spin operators for each qubit a fa, and the raising and lowering pseudo- 
spin operators cr^ and cr^, where each of these operators are given respectively by 



^•)=(|2)(2|- |1)<1|)_, 
2)(1| 



and finally 



(?) 



|1)(2 



(4) 
(5) 

(6) 



This representation is a second quantization of the qubits. Combining Eq. (3) and Eq. (4), the 
j — th qubit Hamiltonian can be written as 



H 



(j) 

D 



*&) + ol"F ) +a>?' ) 



2 v> 2 

where the energy gap between the energy eigenstates of the j — th qubit is given by 



(7) 



(8) 



Shifting the zero of energy to + w^) for each qubit, the j — th qubit Hamiltonian given by 

Eq. (7) can be rewritten as 



H 



(j) 

D 



(9) 



Note that the pseudo-spin operators <r^, and <r^ satisfy the standard angular momentum 
commutation relations corresponding to spin | operators, i.e., 

(10) 
(11) 
(12) 





_2 







and finally 



-2 a 
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To introduce the coupling between the bosonic field and the qubits, let us assume a single quantized 
mode of a bosonic field with a linear coupling with the qubits. Therefore, the Hamiltonian of the 
j — th qubit Hjp, with the contribution of the single quantized mode of the bosonic field Hs, and 
the interaction Hamiltonian Hj^ is given by 

Is <8> H { i ] + H S ® I B + = 

Is ® -y- ct z U) + lu V b <g> I B + g [b + V) <g> + , (13) 

where the second term in the Eq. (13) has the contribution from the single quantized mode of 
the bosonic field Hamiltonian and the last term is the interaction Hamiltonian of the j — th qubit 
with the one-mode quantized field. In the above equation g is a small coupling constant between 
the qubit and the single mode of the bosonic field. The b and fe' are the boson annihilation and 
creation operators of mode excitations that satisfy the usual commutation relation rules. The 
generalization to N qubits is described by 

N N 

I S ®J2 H d +H S ®I B + Y. H i j) = 



3=1 3=1 

JL 

I I " 3 I 



N _ N 



7 = 1 z V^V 1=1 



This Hamiltonian is a possible generalization of the Dicke model. The interaction Hamiltonian 
given by Eq.(14) is simplified if we ignore the counter-rotating terms. Considering this approxi- 
mation we have 

N N 
I s ®Y. H D+Hs®lB + Y. H< i j) = 



3=1 3=1 

(7 (i) -f- UJq 00 Q9 iB i 

3=1 



.7=1 Z ViV =1 



This Hamiltonian describes the original Dicke model in the situation = f2,Vj. This approx- 
imation is known as the rotating-wave approximation. In the rotating-wave approximation we 
ignore energy non-conserving terms in which the emission (absorption) of a quantum of a quan- 
tized field is accompanied by the transition of one qubit from its lower (upper) to its upper (lower) 
state. The rotating-wave approximation ignores terms in which the j — th qubit raising (lowering) 
operators multiplies the mode of the field creation (annihilation) operator. Although the rotating- 
wave approximation has been used to describes the interaction of radiation with matter, as we 
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stressed, there are many situations where we can not assume the above mentioned approximation, 
as for example if we go beyond the dipole approximation, to describe Casimir-Polder forces [48], 
the Lamb-shift and also when atomic systems generate high intensity fields. 

Situations where we can not assume the rotating-wave approximation has been also extensively 
discussed by one of the authors. First Svaiter and Svaiter [49] [50] evaluated the transition rates 
of a two-level atom in different kinematic situations assuming a weak coupling between this small 
system and a real massless scalar field. These authors studied the Unruh-Davies effect [51] [52], 
where an accelerated two-level atom measures a thermal spectrum even if the field is prepared 
in the vacuum state. The origin of this effect is the presence of an event horizon that change 
virtual processes in real ones. In the Ref. [53] still studying the same two-level system coupled to 
a real scalar field, Ford et al assumed the presence of two infinite perfectly reflecting plates which 
change the vacuum fluctuations associated to the quantized bosonic field. The image method 
and the imaginary time formalism [54] [55] [56], was used to study radiative processes at finite 
temperature. In the Refs. [57] [58] one of the authors continue to investigate radiative precesses 
associated to the Unruh-Dewitt detector [59], in interaction with a massless scalar field. Being 
more precise, in the Ref. [58], it was presented the detector's excitation rate when it is uniformly 
rotating around some fixed point, when the scalar field is prepared in the Minkowski vacuum. 

Another model, where the behavior is quite interesting from the mathematical and physical 
point of view is the one where the coupling between the N qubits and the bosonic field in a lossless 
cavity is intensity dependent. We have 



Going back to the case of only one qubit, the two-level atom coupled to a single mode quantized 
electromagnetic field is known as the Jaynes-Cummings model. The Jaynes-Cummings model for 
one qubit, can be written a as 



Practical realization of this model in the laboratory is not a problem. First we consider only one 
two-level atom in the presence of a quantized bosonic field in a cavity. Assuming that the frequency 
c^o of one of the cavity modes is near-resonant with the energy gap Q of the two-level atom, such 
situation generates the following physical model. The two-level atom effectively interacts only 
with that mode, and all the other bosonic modes do not couple with the two-level atom. 



Is ® E Hff + H s ® I B + E Hi = 






(17) 
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This model provides a useful means of studying nonclassical effects in the interaction between 
fields and matter, as for example the phenomenon of collapses and revivals of the Rabi oscillations 
in a quantized field that is not in a pure number state [60] [61] [62] [63]. It is important to point 
out that the behavior of the two- level atom depends if the single field mode is quantized or not. 
In the case of a classical field the model presents atomic population inversion with monotonic 
periodic oscillations. Assuming that the single field mode is quantized the system presents two 
different dynamical behavior. If we prepare the atom in the ground state and the field mode in a 
pure number state, it can be shown that the expectation value of the pseudo-spin operator <r z is 
the same as in the case of a classical field. If we prepare the atom in the excited state and the field 
mode in a pure number state, the atom periodically returns to the upper state with a definite Rabi 
frequency. For the particular case of the atom in an empty cavity (zero occupation number i.e., 
n — 0), the two- level system undergoes periodic, reversible spontaneous decay. Finally, preparing 
the field mode in a coherent state, the dynamics of the two-level system present quantum collapse 
and also revival effects. 

It is possible to generalize the Jaynes-Cummings model, still using the rotating-wave approxi- 
mation in the following way 

I s ® Hff + H S ® I B + H\ 3) = 

Is ® e z (j) +uj bH <g> I B + g(b 2 <g> <j ( + + (tf) 2 <g) (7 ( - } ) . (18) 

In the literature the above model is known as the two-photon Jaynes-Cummings model [64] [65]. 

The same idea that has been used in the Eq. (16) can also be implemented in the Jaynes- 
Cummings model. The Hamiltonian for the Jaynes-Cummings model with a intensity-dependent 
coupling constant where we are still using the rotating-wave approximation is given by [66] [67] 
[68] 

Is <g> H%> + H S ® I B + = 

Is ® a \j) + 6 f b <g> I B + g(b (6 f 6)^ <g> cr+ } + 6 f (tfb)^ <g> cr^) . (19) 

So far we have discussed N qubits interacting with a single quantized mode of the bosonic 
field. Our aim is now to discuss the interaction of a system of N identical qubits with energy gap 
(Q = cu2 — ui), with an infinite number of harmonic oscillators which defines the reservoir. Let a\ 
and ak be the creation and annihilation operators of the k — th harmonic oscillator of frequency uj^. 
These creation and annihilation operators satisfy the standard commutation relations [a' k , a q ] = 8 kq , 

[a k ,a q ] = and also [a|.,at] = 0. The total Hamiltonian, i.e., the Hamiltonian of the combined 
system of the reservoir and the N identical qubits interacting with the reservoir reads 

Q N _ N 

H = Ir® a z U) + J2^ k ala k <g> I s + -^EE^ ® or J) + a I ® %)) ■ ( 20 ) 

Z j=l k Viv j=1 k 
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In the Eq. (20) the first term in the right side is the free Hamiltonian of of N identical qubits, 
the second term is the free harmonic oscillators reservoir Hamiltonian and finally the third term 
is the interaction Hamiltonian between the reservoir and the N identical qubits. Note that we 
shift the zero of energy for each qubits, as we did before, and we are assuming the rotating-wave 
approximation, where g{\fN) ~ x is the j-th qubit, k-th harmonic oscillator coupling constant. 
We also can use a different interaction Hamiltonian to study the influence of decoherence in 
quantum computers as was introduced by Di Vicenzo [5]. This author proposed the following 
model describing a system of one qubit coupled to a reservoir of harmonic oscillators: 

I R ® H S + H R ® I S + H T = 

Ir ®-<x z + E^4 a fc ® I s + gJ2{ a k + a k) ® ° z , (21) 
1 k k 

where Q is the usual energy level spacing of the qubit, a\ and are respectively the boson creation 
and annihilation operators associated to the harmonic oscillators. Note the particular coupling 
between the reservoir and the qubit. It is well known that this model allow for an exact analytic 
solution and also exhibits the destruction of quantum coherence without decay of population. 
There are two straightforward generalizations for this model given by Eq. (21). The first one is 
to introduce N identical qubits and the total Hamiltonian for the composed system reads 

I R ® H S + H R ® Is + H T = 

Q JV N 

!r ® itE 4) + E^4«fc ® J s + ^=EE (4 + «*) ® <AiY ( 22 ) 
z j=l k V iV j=1 k 

Other generalization for the Di Vicenzo model is to introduce a arbitrary mode-dependent coupling 
constant [69] . Therefore we have the following model describing a system of one qubit coupled to 
a reservoir of harmonic oscillators: 

I R <g> H s + H R ® I s + Hi = 

Ir ®7:cr z + E^fc4«fc ® Is + J2{ X k a l + K a k) ® o\ (23) 
1 k k 

Working in this line, the Hamiltonian describing a system of N qubit coupled to a reservoir of 
harmonic oscillators reads [70] 

Ir ® H S + Hr <g> I s + H T = 




where describes the coupling between the j-th qubit with the reservoir. 
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After this briefly discussion of many models describing the two-level atoms-Bose field coupled 
system we would like to go back to Eq.(20). Going back to Eq. (20), another possibility is not as- 
sume the rotating-wave approximation in the interaction Hamiltonian. Without the rotating-wave 
approximation, the interaction Hamiltonian between the N qubits and the reservoir of harmonic 
oscillators reads 

N 

Hi = -j= E E («* + 4) ® + % )) . (25) 

Our aim is to show that the generalized fermion Dicke model where we are not assuming the 
rotating-wave-approximation also present a phase transition and present the spectrum of the 
collective bosonic excitations of the model. In the next section we will study the Hilbert space for 
the fermion Dicke model. The main points we wish to develop in the next section is the functional 
integral for the generalized fermion Dicke model and also discuss the non-analytic behavior of 
thermodynamic quantities at some critical temperature. 



3 The Hilbert space for the fermion Dicke model 

It should be emphasized that there are different routs to proceed in the study of cooperative 
spontaneous emission of a radiation field in the Dicke model. The first one is to introduce collective 
atomic operators 

N 

J' = T,°h> (26) 
i=i 

N 

^ + = E4)' ( 27 ) 
i=i 

and 

N 

J- = T,°(0- ( 28 ) 



i=i 



These operators J z , J + and J~ obey the standard angular momentum commutation relations. 
The Hilbert space of this algebra is spanned by the Dicke states \ j , m) which are eigenstates of J 2 
and J z . This method allow us to describe the N two-level atoms by a single (N + 1) level system. 
A second one method as we stressed before, is to change the atomic pseudo-spin operators of the 
Dicke model by a linear combination of Fermi operators to define the fermion Dicke model. The 
dimensionality of the space where the combinations of Fermi operators act are greater than the 
dimensionality of the space where the pseudo-spin operators act, and therefore we have a problem 
of the elimination of the superfluous states. We start with the Hamiltonian of the Dicke model, 
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H D given by 



N 



N 



H 



D 



i=l y -IV i = i 



(29) 



Note that the Hn Hamiltonian contains the pseudo-spin operators used to obtain a second 
quantized version for the free two-level atoms and also the creation and annihilation operators of 
the single quantized mode of the bosonic field. As we discussed before, the pseudo-spin operators 
a (i) an d CTfii, f° r the same or different qubits obey the standard commutation relations 



'(<) 



a d)> a U) 



and finally 



where 8^2 is the Kronecker delta. 



a (i) d «> 



2a + 



(30) 
(31) 

(32) 



Let us define the Fermi raising and lowering operators a\, oti, fa and fa, that satisfy the anti- 



commutator relations ccjaj + 



Sij and + A T A = Sij. We can also define the following 



bilinear combination of Fermi operators, a\oii — (3}fa, a\fa and finally Pjati. Note that cr*^, afy and 
cr^ obey the same commutation relations as the above presented bilinear combination of Fermi 
operators. This suggests that we can change the pseudo-spin operators of the Dicke model by the 
bilinear combination of Fermi operators 



t, 



and finally 



°fi) — ► ( a l a i - Pi Pi) , 



a< : , — > a\fa , 



'(0 



(33) 
(34) 

(35) 



/From now on we use the usual notation instead of the notation stressing the tensor product space 
of the total Hilbert space of the system. With the substitutions that we defined in the Eq. (33), 
Eq.(34) and Eq.(35) we shall call the resulting Hamiltonian as the fermion Dicke model, i.e., Hp. 
The Hamiltonian of the fermion Dicke model can be written as 



Q N a N 

+ - Pi Pi) + -% £(Wa + 6 f pW 

1 i=i VN i=1 v 



(36) 
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Note that in the Eq. (36) we are still adopting the rotating-wave approximation. We can also 
define the fermion Dicke model, without assume the rotating-wave approximation. The interaction 
Hamiltonian Hi of the single quantized mode of the bosonic field in the presence of the N qubits 
without the rotating-wave approximation is given by 

N 

= ^= E ( &t + h ) («! A + Phi) • (37) 



In effect, this interaction Hamiltonian includes four kind of processes that correspond to the 
absorption or emission of a mode-excitation of the field with the transition of the qubit from its 
lower (upper) to its upper (lower) state. Therefore the interaction Hamiltonian includes virtual 
processes or the vacuum fluctuations contributions. A simple generalization of Eq. (37) that we 
are adopting in the paper is to introduce two different coupling constants, the first one for the 
rotating-wave terms and another one for the counter-rotating terms. We are The advantage of 
this method is that it is possible to identify the contributions of the virtual processes and real 
ones in the formation of the bosonic condensate. 

After this discussion, let us analyze the Hilbert space of the fermion Dicke model. The fermion 
Hilbert space for each atom is four dimensional. The Hilbert space for the i-th atom is generated 
by the following term vectors. The first one is the vacuum state 

00 = 10,0),. (38) 

Next, we can define more two vectors applying the operators a\ and 0\ on the vacuum state, i.e., 
a\(j)Q and 0\<f>o. Therefore we have 

0410,0)* = 1 1,0), (39) 

and 

p\\ 0,0),= | 0,1),. (40) 
Finally, the bilinear combination of the Fermi operators (ot\0\ ) acting on the vacuum state yields 

c4$|0,0), = |l,l),. (41) 

The vectors 1 1, 0), and | 0, l)j generate a two-dimensional subspace, which is characterized by the 
following condition with the bilinear combination of Fermi operators 

(ala i + (3j(3 i )\0,l) i = | 0, 1), (42) 

and 

(<4a i + p}p i )\l,0) i = \l,0) i . (43) 
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We define Ni = (ajccj + 0\fli) as the fermion number operator acting in the Hilbert space 
corresponding to the i-th atom. With the four dimensional Hilbert space we can construct a 
physical Hilbert space and a nonphysical one generated by the vectors | and | respectively. 
Therefore we have 

ci|o,i) i + 4|i,o) i = |*) i 



(44) 



and 



d[\0,0)i + d l 2 \ 1,1), = !$),. (45) 

Next we will obtain a formula which connect the partition function of spin Dicke model to 
the partition function of the fermion Dicke model. Using the definition of the fermion number 
operator of all qubits, and that iVj = (ajctj + p[fli) we have 



N = Ni + ^Nj. 



(46) 



Since the number operator and the Hamiltonian operator commute and using the same notation 
we have 

H F = H i + Y,H j . (47) 

Using the fact that 

fT i |0,0) i = J H i |l,l> i = 0, (48) 

we have that 

Hi\$)i = 0, (49) 

where | $), is the general vector state in the nonphysical subspace of the i-th qubit. The trace 
over the nonphysical states of the i-th qubit vanishes. Therefore we have 



| exp 



i{$ I exp 



$) t = 0. 



Therefore we have 
Tr exp 



(-irTr phys e W [-(3H F \ = {-i)"Tr exp[-(3H a ] , 



(50) 
(51) 

(52) 



and we can present the relation between the partition function of the spin Dicke model to the 
partition function of the fermion Dicke model. A simple formula is given by 



Tr exp [-PH a \ = i N Tr exp 



( nrr ™ 



N 



(53) 
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According the Eq. (53), to study the spin Dicke model we can use the fermion Hamiltonian Hp 
adding i.e. the fermion number operator with pure imaginary coefficient (chemical potential). 
All the standard diagrammatic technique for the Fermi system can be generated using the Fourier 
representation for the Green function given by 

G°M = 1 —— = - , (54) 

iujf — e + fi iujf — e — || 

where ojf = + 1/2) is the fermion Matsubara frequency. After the above discussion we can 
consider the problem of define the partition function of the fermion Dicke model defined by Zp. 



4 The functional integral for the generalized fermion Dicke 
model 

After the above discussion we can consider the problem of define the partition function Zp of 
the generalized fermion Dicke model. First let us define the Euclidean action S of this model, 
which describes a single quantized mode of the field and the ensemble of N identical qubits. The 
Euclidean action S is given by 

S = J Q rfr(r(r)-6(r)+^(«*(r)-^(r)+A*(r)-A(r))) - jf drH F (r) (55) 
where Hp is the full Hamiltonian for the generalized fermion Dicke model given by 
H F = u b*(r)b(r) + 5 f; («;(r) a! (r) - /3*(r)/3,(r)) + 
+ j= E(aKr)p i (r)b(r) + a< (r) (3*{r) b\r) ) + 

n N / \ 
+ E («<W^ r W^) + «*(r)^(r)6*(r)j. (56) 

Note that we are introducing two coupling constants, g\ and g<i, for the rotating and anti-rotating 
wave terms respectively. As we discussed before, the main reason for this is that we are interested 
to identify the contribution of the real and virtual processes in the phase transition with the 
formation of the condensate. Let us define the formal quotient of two functional integrals, i.e., the 
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partition function of the generalized fermion Dicke model and the partition function of the free 
fermion Dicke model. Therefore we are interested in calculate the following quantity 



(57) 



Z Fo Jld V ]eSo> 



where S = S(b,b*, a, a* , (3, /3') is the Euclidean action of the generalized fermion Dicke model 
given by Eq. (55), Sq = So(b, b*, a, cr , f3, /r ) is the free Euclidean action for the free single bosonic 
mode and the free qubits, i.e., the expression of the complete action S taking g± — g 2 — 
and finally [drj\ is the functional measure. The functional integrals involved in Eq. (57), are 
functional integrals with respect to the complex functions b*(r) and b(r) and Grassmann Fermi 
fields a*(r), «i(r), 13* (t) and /^(r). Since we are using thermal equilibrium boundary conditions, 
in the imaginary time formalism, the integration variables in Eq. (57) obey periodic boundary 
conditions for the Bose field, i.e., b{(5) = 6(0) and anti-periodic boundary conditions for Fermi 
fields i.e., ccj(/3) = — ccj(0) and Pi((3) = — A(0). 

The free action for the single mode bosonic field S (b) is given by 



Then we can write the action S of the generalized fermion Dicke model, given by Eq. (55), using 
the free action for the single mode bosonic field So(b) given by Eq. (58), plus an additional term 
that can be expressed in a matrix form. Therefore the total action S can be written as 




(58) 




i=l 



(59) 



where Piij) is a column matrix given in terms of fermion field operators 




(60) 



and the matrix M{b* , b) is given by 




(61) 
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These fields b(r), a^r) and A(r) can be written as a Fourier expansion. Therefore we have 

^) = r 1/2 EK^ r , (62) 

U! 

and 

Pi(T)=p- 1/2 EPi(P> i * T - (63) 
p 

Since the field 6(r) obeys periodic boundary conditions, and the fields a^r) and /^(t) obey anti- 
periodic boundary conditions, we have that uj = ^ and p = ( 2n + 1 ) 7r ; where they are the boson 
and fermion Matsubara frequencies respectively. Substituting the Fourier expansions in the free 
action given by Eq. (58) we get 

S (&)=5>^-wo)&*H&H. (64) 
Using the same procedure in the Eq. (59) allow us to write 

S = S (b) +EE pKp) M P,(b*, b) Pi (q) , (65) 
p,q i=i 



where the matrix M pq (b*, b) is given by 
/ 

M p(? (&*,5) = 



V 



(ip + fi/2)5 pg (iV/3)-V2 ^ _ p ) + g2 6(p _ g) 

(iV^)- 1 ^ ^ & (p _ 5 ) + ^ _ p )) ( ip _ n/2)5 pq 



(66) 

Using the above results, the ratio between the two functional integrals Z and Z , i.e., ^- is given 
by 

f[dr,(b)] exp faiu; - uj )b*(uj)b(uj)) f [d V (p)\ exp(]T £ pl(p) M pq (b*, b) Pi {q) 

^ w ' J \p,q i=l 



/[di;(6)] expfaiu; - J [dn{p)\ exp(£ £ pj(p) M pg (0, 0) p,(g) 

^ w ' ^p,q i=l 



(67) 



where the functional measures and [drj(p)] in the above equation are defined respectively 

by 

[dri(b)]=l[db(u)db*(u) (68) 
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and 

[dv(p)] = Ud Pi (p)dp\(p). (69) 

i,p 

We need impose cutoffs over the boson and fermion Matsubara frequencies on these measures. 
This procedure is necessary to be sure that the ratio between the two functional integrals given 
by Jr- does not diverge. After all , at the end, we must take these cutoffs to infinity. In order to 
define the effective action associated to the bosonic mode, we integrate out the fermionic degrees 
of freedom. The integrals with respect to the Fermi fields are Gaussian and we may integrate over 
these Grassmann variables. This procedure yields 

/ [dri(p)] exp(£ J2 pKp) M Pq (b*, b) Pl {q)) = det N M(b*, b) , (70) 

J \p,q i=l > 

where the matrix M is a block matrix of the following form 

where I is the identity matrix and the components of matrix Q are 

Q Pq = 9ib(p-q)+g2b*(q-p). (72) 

The following change of coordinates can simplify our calculations. Let us change variables in the 
following way: 

1/2 

(73) 

and 




&». (74) 



We must note that Eq. (74) is not conjugate of Eq. (73). Nevertheless, is not difficult to 
justify this transformation, if we introduce polar coordinates instead of b(u>), b*(u) : b(u) = 
{p{u)) 1 / 2 e^' J} \ b*(u) = (p(c<j)) 1 / 2 e" i<?f '^ and then perform a complex rotation of the integration 
counter when integrating with respect to p(uj) : p(uj) — > p{uj)[n / (u — iuj)] 1 / 2 . Its easy to see that 
after these changes of variables the denominator of the Eq. (67), turns out to be equal to unity 

J [d77(6)]exp(-7r£r(u;)&H) = 1, (75) 
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so we can express the ratio by the integral 



Z 



=r = J[dT,(b)]exp(s ef ,(b)}, (76) 

where S e ff(b) is the effective action of the bosonic mode which is given by 

S eff = -7T b*{uj)b{uj) + N\ndet(I + A). (77) 

The determinant in the above equation is given by 

det(I + A) = det(^M- 1/2 (0,0)M(6*,6)M- 1/2 (0,0)) (78) 
and the matrix A is defined as follows 

•4 - [ ;, : ! . (79) 




In the equation above the quantities B and C are matrices with components given by 

~ v^v) r + 2 j 1 ^ _, + I 19 - 1 j (80) 



and 



/37V/ v 2/ V - i(p - q) yjuo - i(q-p) 



In Eq. (76) we may go to the limit wg, oup — > oo and the instead of a formal quotient of two 
infinite functional integrals we shall have only one finite functional integral. This representation 
turns out to be very useful for obtaining the asymptotic formula for Z/Z at large N. There exists 
only one stationary phase point at > If f3~ l < f3~ l , we have a circle of a stationary phase 
| 6(0)| 2 = po? b(u>) = b*(u>) = 0, if u ^ 0. There also exists an interpolation formula between these 
asymptotes. The presence of degenerate vacua is a feature of states with spontaneous symmetry 
breaking. As we will see, gapless excitation will appear. 

We shall investigate the integral given by Eq. (76) for temperatures that satisfy > 
First of all let us show that this integral converges. We use the estimate 

| det(J + A)\< exp (Re (trA) + ^ tr(AA^ . (82) 



18 



where Re (trA) means the real part of trA. The matrix A have the form given by Eq. (79). 
Therefore we find that trA = and tr(AA*) = tr(BB^) + tr(CC f ). Therefore, we obtain the 
estimate 

< f[dr](b)} exp( -TrJ^b* (co)b(cu) + Ntr(BB^) + Ntr(CC^)) , 

< J [d V (b)\ exp^-vr 6» (l - oqH) b(u) + 

+ (b(u) c oH K - to) + 6» c (lo) b*(- u) ) V (83) 



where the a (u) and c (cj) are given respectively by 

g\ + g\ 



and 



a (u}) = 



E 



i 



i 



(3{UI + W ')V2 ^ (Hi + ,2)1/2 + p 2)l/2' 

^0 9l92 I ' 



(84) 



(85) 



/?K 2 + ^) p _v = . (¥ + 1 2 ) 1 ' 2 (¥ + p 2 ) 1/2 ' 

Using the measure given in Eq. (68), we have that < F, where F = F 1 F 2 and where Fi and 
i<2 are given by 



Fi = y db{o)db*(o) 

tt6*(0) (l - a (0))6(0) + 7T (6(0) c (0) 6(0) + 6*(0) c (0) 6*(0) 



cxp 



(86) 



and 



F 2 = I [] db{uj)db*{uj)db{-uj)db*{-uj) 
->o 

2 b*(u) (l - a H ) b(u) - vr £ 6*(-w) (l - a„M ) 6(-w) + 



cxp 



lu>0 



uj>0 



2tt (b(ou)c (u)b(-u) + b*{uj)c (uj)b*{-uj)) 

oj>0 ^ ' 



(87) 



Note that in the case of the generalized fermion Dicke model we obtained a Gaussian integral that 
mixtures positive with negative frequencies. A straightforward calculation gives that the ratio 
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obeys the following inequality 

1 - oo(0) + 2c (0)) (l - ao(0) - 2c (0) 



Z 

— < 



■1/2 



n 

w>0 



^1 - a (u;) + 2c (cj)^ ^1 - a (a;) - 2c (^) 



- i 



(88) 



In a similar way like Popov and Fedotov [20] proved, for the case of rotating wave approximation, 
we have that, < a (uj) + 2c (cj) < a (0) + 2c (0) and a (0) + 2c (0) = 0(cj~ 2 ln a;). Therefore 
if a (0) + 2c (0) < 1, then Eq. (88) guarantees convergence of the expression J^. The condition 
a (0) + 2c (0) = 1 is the equation for the transition temperature, then we have 



ao(0) + 2c (0) = 



.91 + 92 



-(¥ 



= 1 



The inverse of the critical temperature /3 C is given by 

4 , _-, / f2a>o 



& = - tanh 1 (- p-r^) ■ 

" V(^l+^2) 2 / 



(89) 



(90) 



Note that there is a quantum phase transition where the coupling constants g\ and g 2 satisfy 
9i + #2 — (^o^) 5 - F° r larger values for (gi + (? 2 ) the system enters in a superradiant phase. 

To calculate the asymptotic behavior of the functional integrals at temperatures that satisfy 
> [3" 1 , we can do the following approximation 



det N {I + A) = det N (I + BC) -> exp(iVtr(5C)) 



(91) 



This substitute can be done and we can estimate the error if we divide all the functional space 
into two domains C\ and C 2 



Denoting 



tr[(BC)(BCy] < (AN)- 1 i-> d , 
tr[{BC){BCy] > {AN)- 1 i-> C 2 . 

= det ^(J + A) - exp(iVtr(flC)) , 



(92) 
(93) 

(94) 
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for the ratio -J^, we have the following identity 



Z 
~Zi 



- = [ [dr]{b)}e^(-nJ2b*(u;)b(uj) + Ntr(BC)) + 
o J V w / 

+ f c [d V (b)]K N exp(-irJ2b*(u;)b(u;f) + 
+ J c [dr]{b)}K N exp(-7r^2b*(uj)b(u 



(95) 

The first integral of the above equation is Gaussian, let us define it by I . We use the Eq. (80) and 
the Eq. (81) in order to calculate the trace of BC, i.e., tr(BC) which is present in the expression 
Iq. A simple calculation gives 

Jo = J [dri(b)\ exp^-vr ]T &» (l - a(w)) b(u>) + 
+ 7i J2 c(w) &(-<<;) + 6*(w) c(w) &*(- cc 



(96) 



where a (a;) and c(u;) of above equation are given respectively by 

(g\ (fi - iu) - l + 5-| (fi + iw) _1 



a (a;) 



and 



c(w) = 



(u - iu;) 



tanh 



(co 2 + CJ 2 ) 1 /2(^2 + U 2y 



tanh 



pn 



(97) 



(98) 



Note that to recover the result obtained by Popov and Fedotov [19] [20] we have only to assume 
g 2 = 0. In this case, we have that c(u>) = 0, which simplify the integration over the mode of the 
bosonic field in Eq. (96). Making the integration we obtain that I is given by 



/ =n(i-a 



(99) 



After this observation lets go back to the general case, where g\ and g 2 take arbitrary values. The 
expression I given in Eq. (96) is a Gaussian integral, this expression is similar to the integral 
given in Eq. (83), so following the same steps we get that 



j = i (u = o) n 

u;>0 



c(uj) 2 - (l - a(u)^ (l - a(-w)) 



(100) 
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where Iq(u = 0) is the contribution of the condensate given by 



h{oo = 0) 



1 - a(0) + 2c(0)) (l - a(0) - 2c(0)) 



1 -1/2 



(101) 



It is possible to estimate the error of I which is given by the two last terms of Eq. (95). For 
details see the reference [8]. The errors depend on order iV _1 . Therefore can be written as 



Z 



1 - a(0) + 2c(0) 1 - a(0) - 2c(0) 



n 

w>0 L 



1 - a(w)) (l - a(-w)) - c 2 (cj) 



-1/2 



+ 



(102) 



+C(7V- 1 ), 

therefore in the limit (N — > oo) the equality = I is a good approximation. To find the 



collective excitation spectrum we have to use the equation 

c 2 (oj) - (l - (l - a(-w)) =0, 

and making the analytic continuation —>£'), we obtain the following equation 



(103) 



1 = - 



4 i 4 

91 + 91 



K - £ 2 )(ft 2 - E*) 
( 1 



tanh 2 ' 



+ 



9\9l 



+ 



4ft' 



(a; 2 — .E 2 ) \ (ft — E) 2 (ft + £) 2 (ft 2 -£ 2 )' 



+ 



4 2 (ft - £ ) _1 + «? 2 2 (ft + E)~ l + fl 2 (ft + i?)- 1 + £ 2 (ft - E) 



-1 



(c^o - E) 



(u> + E) 



+ 



tanh 



V 4 



(104) 



Solving the above equation for the case f3 1 = (3 C 1 , we find the following roots 

E 1 = (105) 

and 

E2 = f <M» + ".)' + (n-".)' ^ (106) 



(#1 + #2) 
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Its low energy state of excitation is a Goldstone mode. Now, let us present the critical temper- 
ature and the spectrum of the collective bosonic excitations of the model with the rotating-wave 
approximation, where gi ^ and g 2 = 0. The result obtained by Popov and Fedotov is recovered, 
where the equation 

a(0) = 1 (107) 

and 

4-tanhf^Ul, (108) 



LU tt V 4 

give the inverse of the critical temperature, f3 c . It is given by 

) 9 c = itanh- 1 (^). (109) 

In this case, also there is a quantum phase transition, i.e., a zero temperature phase transition 
when g 1 = (cu Q) "2 . It is interesting to point out that there are two different ways to analyze 
the phase transition. The first one is to follow the non-analytic behavior of the thermodynamic 
quantities as a function of temperature. A different way is to follow the non-analytic behavior of 
the thermodynamic quantities as a function of the coupling constant strength. Working in this 
second approach, we may expect that for large coupling constant g± there is a superradiant phase. 
The spectrum of the collective Bose excitations in this case is 

E 1= 0, (110) 

and 

E 2 = n + u . (111) 

Now we will show that it is possible to have a condensate with superradiance in a system of N 
qubits coupled with one mode of a Bose field where only virtual processes contribute. In the pure 
counter-rotating wave case, i.e., g\ = and g 2 7^ 0, the inverse of the critical temperature, f3 c is 
given by 

and the spectrum of the colective Bose excitations given by 

E 1= 0, (113) 

and 

E 2 = \Q-u \. (114) 

A comment is in order concerning the Bose excitations spectrum. In the both cases: using or 
not the rotating-wave approximation, there is a phase transition. In the case of the rotating-wave 
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approximation g\ 7^ and g 2 = 0, there is a Goldstone mode (E = 0). In the pure counter- 
rotating wave case g 1 = and g 2 7^ 0, also there is a Goldstone (gapless) mode. The existence 
of Goldstone modes and the energy of the other mode was presented for both above mentioned 
cases. The spectrum in the general case is given by the Goldstone mode and also by a non-zero 
energy mode given by Eq. (106). It is interesting to stress that we obtained a critical behavior in 
both situations (#1 7^ 0, g 2 = and g 1 = 0, g 2 7^ 0), where the condensate has Goldstone (gapless) 
modes, with a superradiant state. Therefore we show that it is possible to have a condensate with 
superradiance in a system of N qubits coupled with one mode of a Bose field where only virtual 
processes contribute. 

We would like to point out to the reader that there is an analogy between the generalized Dicke 
model and the Landau-Ginzburg model in the theory of phase transition. To describe the second- 
order phase transition Landau introduced the concept of order parameter. This concept enabled 
Landau to write the free energy of any system in the form of an expansion in powers of the order 
parameter. This expansion characterizes order in the system as a whole. A generalization of this 
model accounts for the local space variation, i.e., to take into account possible local fluctuation 
of the order parameter. In the Landau theory we ignore fluctuations as the same way that in the 
rotating-wave approximation in the Dicke model. Going along this line we can make a parallel 
between the Landau-Ginzburg model and the generalized Dicke model. In concluding we remark 
that the model with pure virtual processes is enable to generate a condensate. 



5 Conclusions 

In the present paper we are first discussing Hamiltonians describing a quantized bosonic field 
interacting with two-level atoms. Without assuming the rotating-wave- approximation, with the 
coupling constants g\ and g 2 for rotating and counter-rotating terms respectively, we define the 
generalized fermion Dicke model. Changing the atomic pseudo-spin operators of the model by a 
linear combination of Fermi operators we define the fermion generalized Dicke model. Studying 
the case where identical two- level atoms act as a thermal reservoir (N — > 00), we investigate the 
thermodynamic of the generalized Dicke model using the path integral approach with functional 
integration method. We are considering the question of how does the counter-rotating terms of 
the interaction Hamiltonian contribute in the formation of the condensate with a superradiant 
phase transition in the model. 

We study the nonanalytic behavior of thermodynamic quantities of the generalized model for 
both situations, evaluating the critical transition temperature and presenting the spectrum of the 
colective bosonic excitations, for the case gi ^ and g 2 — 0, g\ — and g 2 7^ and also in 
the general case. Our result show that it is possible to have a condensate with superradiance 
in a system of N qubits coupled with one mode of a bosonic field where only virtual processes 
contribute. It is important to realize that the energy of the non-Goldstone mode in Eq.(lll) is 
always larger than the energy of the non-Goldstone mode of Eq. (114), i.e., in the system where 
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the condensate appears due to the virtual processes. Our conclusion from the above results is that 
both processes, real and virtual ones give different contributions to generate the condensate. 

An important question is the way of practical realization of the generalized Dicke model in the 
laboratory. As was stressed by Dimer et al [71] it remains as a challenge to provide a physical 
system where the counter-rotating terms are dominant. Those authors proposed that in cavities 
with the N qubits, only one mode of quantized field and classical fields (lasers), it is possible 
to obtain a physical system that corresponds to the generalized Dicke model. Also, has been 
discussed in the literature the possibility of control the relative importance of the counter-rotating 
terms in the Jaynes-Cummings model in the laboratory using a ion trap [72]. Another mechanism 
to explore the importance of virtual processes was proposed by Ford [73] and Ford and Svaiter 
[74] [75], where the possibility of amplification of the vacuum fluctuations was discussed. These 
authors studied the renormalized vacuum fluctuations associated with a scalar and electromagnetic 
field near the focus of a parabolic mirror. Using the geometric optics approximation these authors 
found that the mirror geometry can produce large vacuum fluctuations near the focus. 

It is not our proposal to discuss the practical realization of the generalized Dicke model in the 
laboratory, but only analyze the thermodynamic of the model using the path integral approach 
with functional method. An evidence in favor of our results is an experiment where it is possible 
to control the importance of the counter-rotating terms in the generalized model in such a way 
that an ideal g± ~ situation is achieved. Experimental observation of the superradiant phase 
in this situation will improve our understanding of this phenomenon, when the atoms radiates 
spontaneously, at a radiation rate much higher than would be expected from an ensemble of 
independently radiating atoms. 

There are many different continuations for this paper. The first one is to investigate the model 
introduced by DiVicenzo, defined by Eq. (22), at finite temperature, using also functional integral 
methods. Also the model defined by Eq. (16), with a coupling between the N qubits and the 
single quantized mode of a bosonic field, which characterize a intensity-dependent coupling can be 
discussed using the conventional technique of the path integral and functional integration method. 
A similar calculation presented in this paper could be carried out for this model, although the 
calculations would be somewhat more laborious that the presented in the paper. The generalization 
of this model with the introduction of the counter-rotating terms deserves further investigations. 

An outstanding question is the presence of the quantum phase transition and also the chaotic 
behavior for a finite N, where a crossover between Poisson and Wigner-Dyson behavior in the 
energy level-spacing statistics appears. The question that has been discussed in the literature 
is the distribution function of spacing of adjacent eigenenergies in a large number of systems. 
For systems where the classical dynamics is integrable, a Poissonian nearest neighbor spacing 
distribution is expected, i.e., the spacing rule for random levels. Some systems, the eigenenergies 
do not follow the Poisson law, but behave as the eigenvalues of a random matrix taken from 
a suitable ensemble. The set of all real random matrices with matrix elements obeying some 
distribution function defines the Gaussian orthogonal ensemble. Other ensembles are the Gaussian 
unitary and the Gaussian sympletic ensemble respectively. The chaotic behavior for a finite N in 
the above discussed model is under investigation by the authors. 



25 



6 Acknowlegements 

We would like to thanks Alexis Hernandez for useful discussions. This paper was supported by 
Conselho Nacional de Desenvolvimento Cientifico e Tecnologico do Brazil (CNPq). 

References 

[1] R. H. Dicke, Phys. Rev. 93, 99 (1954). 

[2] D. Meshede, Phys. Rep. 211, 201 (1991). 

[3] P. Meystre, Phys. Rep. 219, 243 (1992). 

[4] E. T. Jaynes and F. W. Cummings, Proc. Inst. Elect. Eng. 51, 89 (1963). 

[5] D. P. DiVicenzo, Phys. Rev. A51, 1015 (1995). 

[6] H. Dekker, Phys. Rep. 80, 1 (1981). 

[7] K. Blum, Density Matrix Theory and Applications, Plenum Press, N. Y. (1981). 

[8] V. N. Popov, Functional Integrals and Collective Excitations, Cambridge University Press 
(1987). 

[9] A. Zardecki, Phys. Rev. A22, 1664 (1980). 
[10] M. Hillery and M. S. Zubaire, Phys. Rev. A26, 451 (1982). 
[11] E. Schrodinger, Naturwiss 14, 664 (1926). 
[12] R. Glauber, Phys. Rev. Lett 10, 84 (1963). 
[13] R. Glauber, Phys. Rev. 130, 2529 (1963). 
[14] R. Glauber, Phys. Rev. 131, 2766 (1963). 
[15] E. C. G. Sudarshan, Phys. Rev. Lett. 10, 227 (1963). 

[16] P. L. Knight and L. Allen, "Concepts in Quantum Optics", Pergamon Press Inc., N.Y. (1983), 
pp. 54-69. 

[17] J. Kondo, Physica 125B, 279 (1984). 

[18] L. D. Chiang and S. Chakravarty, Phys. Rev. B, 154 (1985). 

26 



[19] V. N. Popov and S. A. Fedotov, Theor. Math. Phys. 51, 363 (1982). 

[20] V. N. Popov and S. A. Fedotov, Sov. Phys. Jetp. 67, 535 (1988). 

[21] J. R. Ackerhart, P. L. Knight and J. H. Eberly, Phys. Rev. Lett 30, 456 (1976). 

[22] P. W. Milonni, J. R. Ackerhart and W. A. Smith, Phys. Rev. Lett 31, 958 (1976). 

[23] P. W. Milonni and W. A. Smith, Phys. Rev. All, 814 (1976). 

[24] B. Fain, N. Cimento B68, 73 (1982). 

[25] J. R. Ackerhart and J. H. Eberly, Phys. Rev. D10, 3350 (1974). 

[26] J. Dalibard, J. Dupont Roc and C. Cohen Tannoudji, J. Phys. 43, 1617 (1982). 

[27] J. Dalibard, J. Dupont Roc and C. Cohen Tannoudji, J. Phys. 45, 637 (1984). 

[28] P. W. Milonni, "The Quantum Vacuum", Academic Press Inc., San Diego (1994), pp. 137-142. 

[29] S. Sachdev, "Quantum Phase Transitions", Cambridge University Press, Cambridge (1999). 

[30] V. B. Kir'yanov and V. S. Yarunin, Theor. Math. Phys. 43, 340 (1980). 

[31] K. Hepp and E. H. Lieb, Ann. Phys. 76, 360 (1973). 

[32] K. Hepp and E. H. Lieb, Phys. Rev. A8, 2517 (1973). 

[33] Y. K. Wang and F. T. Hioe, Phys. Rev. A7, 831 (1973). 

[34] F. T. Hioe, Phys. Rev. A8, 1440 (1973). 

[35] G. Comer Duncan, Rev. A8, 418 (1974). 

[36] T. Holstein and H. Primakoff, Phys. Rev. 58, 1098 (1949). 

[37] L. R. Mecdand and N. Papanicolau, Phys. Rev. B28, 1633 (1983). 

[38] M. Hillery and L. D. Mlodinow, Phys. Rev. A31, 797 (1985). 

[39] C. Emary and T. Brandes, Phys. Rev. Lett 90, 044101 (2003). 

[40] C. Emary and T. Brandes, Phys. Rev. E67, 066203 (2003). 

[41] R. Graham and M. Hoherbach, Phys. Lett. 101A, 61 (1984). 

[42] R. Graham and M. Hoherbach, Phys. Rev. Lett. 57, 1378 (1986). 

27 



[43] C. H. Levenkopf, M. C. Nemes, V. Marvulle, M. P. Pato and W. F. Wreszinski, Phys. Lett. 
A155, 113 (1991). 

[44] L. Muller, J. Stolze , H Leschke and P. Nagel, Phys. Rev. A44, 1022 (1991). 

[45] G. A. Finney and J. Gra-Banacloche, Phys. Rev. A50, 2040 (1994). 

[46] P. W. Milonni, J. R. Ackerhalt and H. W. Galbraith, Phys. Rev. Lett. 50, 966 (1983). 

[47] R. Kullock and N. F. Svaiter, "A Continuously Observed Two-level System Interacting with 
a Vacuum Field", quant-ph/0703167 (2007). 

[48] H. B. G. Casimir and P. Polder, Phys. Rev. 73, 360 (1948). 

[49] B. F. Svaiter and N. F. Svaiter, Phys. Rev. D46, 5267 (1992). 

[50] B. F. Svaiter and N. F. Svaiter, Phys. Rev. D47, 4802 (1993) (erratum). 

[51] P. C. W. Davies, J. Phys. A8, 365 (1975). 

[52] W. G. Unruh, Phys. Rev. D14, 870 (1976). 

[53] L. H. Ford, N. F. Svaiter and M. L. Lyra, Phys. Rev. A49, 1378 (1994). 

[54] R. Kubo, J. Phys. Soc. Jap. 12, 570 (1959). 

[55] P. C. Martin and J. Schwinger, Phys. Rev. 115, 1342 (1959). 

[56] S. A. Fulling and S. N. M. Ruijsennars, Phys. Rep. 152, 135 (1987). 

[57] V. A. De Lorenci, R. De Paola and N. F. Svaiter, Class. Quant. Grav. 20, 4241 (2000). 

[58] R. De Paola and N. F. Svaiter, Class. Quant. Grav. 18, 1799 (2001). 

[59] B. S. Dewitt, Phys. Rep. 196, 295 (1976). 

[60] J. H. Eberly, N. B. Narozhny and J. J. Sanches-Mongradon, Phys. Rev. Lett. 44, 1323 (1980). 
[61] N. B. Narozhny and J. J. Sanches-Mongradon and J. H. Eberly, Phys. Rev. A23, 236 (1981). 
[62] H-I Yoo, J. J. Sanches-Mongradon and J. H. Eberly, J. Phys. A14, 1383 (1981). 
[63] Y. I. Yoo and J. H. Eberly, Phys. Rep. 118, 239 (1985). 
[64] A. Joshi and R. R. Puri, Phys. Rev. A42, 4336 (1990). 



28 



[65] A. H. Toor and M. S. Zubairy, Phys. Rev. A45, 4951 (1992). 
[66] B. Buck and C. V. Sukumar, Phys. Lett. 81A, 132 (1981). 
[67] B. Buck and C. V. Sukumar, J. Phys. A17, 885 (1984). 
[68] V. Buzek, Phys. Rev. A39, 3196 (1989). 
[69] F. Benatti and R. Floreanini, quant- ph/507271. 

[70] G. M. Parma, K,-A. Seiominen and A. K. Ekert, Proc. R. Soc. Lond. A452, 567 (1996). 
[71] F. Dimer, B. Estienne, A. S. Parkins and H. J. Carmichael, Phys. Rev. A75, 013804 (2007). 
[72] J. C. Cirac, A. S. Parkins, R. Blat and P. Zoller, Phys. Rev. Lett. 570, 556 (1993). 
[73] L. H. Ford, Phys. Rev. A58, 4279 (1990). 

[74] L. H. Ford and N. F. Svaiter, Phys. Rev. A62, 062105-1 (2000). 
[75] L. H. Ford and N. F. Svaiter, Phys. Rev. A66, 062106-1 (2002). 



29 



